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We propose methods to generate and manipulate vortex dipoles in a Bose-Einstein condensate 
using Gaussian beams of red or blue-detuned laser. Velocity-controlled vortex dipoles are shown 
to be created and launched by a red-detuned beam and by two blue-detuned beams. Critical 
velocities for the vortex nucleation are investigated. The launched vortex dipoles can be trapped, 
curved, accelerated, and decelerated using Gaussian beams. Collisions between vortex dipoles are 
demonstrated. 
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I. INTRODUCTION 



A vortex dipole is a pair of vortices with opposite cir- 
culations, which is a fundamental excitation in a two di- 
mensional (2D) fluid and plays a pivotal role in thermo- 
dynamic behavior of 2D systems [H, H|. In superfluids, 
vortices are topological defects of the macroscopic wave 
function and a vortex dipole consists of a quantized vor- 
tex and antivortex, whose circulations are quantized to 
±h/m with h and m being Planck's constant and the 
mass of a particle [3|, |4[ . Such a topological structure of 
a quantized vortex dipole has been observed in a Bose- 
Einstein condensate (BEC) of an atomic gas [5[ and in 
an exciton-polariton condensate [6[. Recently, dynamics 
of vortex dipoles have directly been observed in oblate 
BECs 00- 

In the experiments reported in Refs [E|, 0], vortex 
dipoles were created by stirring a BEC with a blue- 
detuned laser beam. When the velocity of the beam 
exceeds a critical velocity, vortex- antivortex pairs are re- 
leased from the density dip produced by the laser beam. 
The created vortex dipoles then move through the BEC 
and exhibit various interesting dynamics |7l-[lQj, which 
can be observed in real time using the technique devel- 
oped by the Amherst group [8|, [9| . Theoretical study on 
a moving potential in superfluids has so far been focused 
on drag for ce Jlll -[l4{ , vortex formation and shedding dy- 
namics [l5l-ll7(, critical velocity [HI, [l9| , scaling laws [20| , 
and potential oscillation [2l|. In addition to these stud- 
ies, it is important to establish methods for controlled 
generation and manipulation of vortex dipoles, which in- 
creases the range of experiments and leads to precise un- 
derstanding of dynamical properties of vortex dipoles. 

The aim of the present paper is to propose methods 
to generate and manipulate quantized vortex dipoles in 
a BEC using not only a single blue-detuned laser beam 
but also a red-detuned beam and multiple beams. We 
will show that dynamics of vortex dipole creation for a 
red-detuned beam (attractive potential for atoms) are 
quite different from those for a blue-detuned beam (re- 



pulsive potential for atoms). We propose methods to cre- 
ate and launch vortex dipoles with a controlled velocity 
and angle using a red-detuned beam or two blue-detuned 
beams. We also show that two red-detuned beams can 
trap a vortex dipole between them. The trajectories of 
the launched vortex dipoles can be curved by both blue- 
and red-detuned beams. The launched vortex dipoles can 
also be accelerated and decelerated by Gaussian beams 
with time-dependent intensity. Using two vortex-dipole 
launchers, we demonstrate collisions of vortex dipoles. 

The present paper is organized as follows. Section [TT] 
provides the formulation and numerical method. Sec- 
tion Hn] studies vortex dipole generation by a red-detuned 
beam and shows that it can launch a vortex dipole in a 
controlled manner. Section HVl investigates dynamics for 
two laser beams. Section [V] proposes methods to ma- 
nipulate launched vortex dipoles and demonstrate their 
collisions. Section fVTI gives conclusions to this study. 



II. FORMULATION OF THE PROBLEM 

We consider a zero-temperature BEC in the mean-field 
approximation. The dynamics of the macroscopic wave 
function ^ is described by the Gross-Pitaevskii (GP) 
equation given by 



2m 



(1) 

where m is the atomic mass, V and V z are the exter- 
nal potentials, and a is the s- wave scattering length. We 
assume that the system is tightly confined in the z di- 
rection and the confinement energy is much larger than 
other energy scales of the system. The wave function 
is thus frozen into ^/(x^y^z^t) = ip{x^y^t)e lEot ^ n ipQ{z)^ 
where ^o(z) and Eq are the ground state and its energy 
for V z (z). Multiplying ^o(z) to Eq. (pQ) and integrating 
it with respect to z, we reduce the 3D GP equation to 
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2D as 



ih 



dt 
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VliP + V(x,y,t)iP + g\iP\ 2 xP, (2) 



where is the 
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2D Laplacian and g = 
is the effective interaction pa- 
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rameter in 2D space. We normalize the 2D wave 

~ —1/2 

function as ijj = n i/j, where no is the 2D atomic 
density in the absence of V. Normalizing the length 
and time by £ = %/ '(mgno) 1 ^ 2 and r = ^/v s = tt/(grio), 
where v s = (gno/m) 1 ! 2 is the sound velocity, Eq. ([2]) is 
rewritten as 



1 dt 



(3) 



where the tilde represents normalized quantity. From this 
normalization, we find that the scaling law holds for the 
mean-field approximation. For example, the dynamics 
for the interaction strength gno with a Gaussian potential 

V = Vo exp{ — [(x — vt) 2 + y 2 ]/d 2 } and that for agn with 

V = aV exp{-[(x - ^vt) 2 + y 2 )/(d/^a) 2 } obey the 
same equation (j3]) for an arbitrary constant a > 0. 

In the numerical simulations, we first prepare the 
ground state for V(x,y,t = 0) by the imaginary time 
propagation of Eq. (|3]), i.e., i is replaced by —1 on the 
left-hand side of Eq. ((3]). We then add a small white 
noise to the ground state to break the exact numerical 
symmetry. Starting from this initial state, we obtain the 
time evolution by solving Eq. (j3j). The imaginary and 
real time propagation is calculated by the pseudospectral 
method, and therefore the periodic boundary condition 
is imposed. We take a large enough space in the simula- 
tions so that the boundary does not affect the dynamics. 



III. VORTEX GENERATION BY A 
RED-DETUNED GAUSSIAN BEAM 

A. Flow around a moving potential 

First, we consider a single Gaussian laser beam moving 
through the condensate, which produces a potential as 



V(x, y, t) = V exp 



(x + vt) 2 + y 2 



d 2 



(4) 



where d is the 1/e width of the laser beam. We assume 
that the laser beam moves in the — x direction with the 
velocity v. The constant Vo is proportional to the inten- 
sity of the laser beam, which is positive for blue detuning 
and negative for red detuning with respect to the reso- 
nant frequency of the electric dipole transition of atoms. 
The detuning is large and the heating of atoms by the 
spontaneous emission is negligible. 

Frisch et al showed that quantized vortex dipoles 
are created around a moving rigid disk (V = +oo inside 
a disk) above a critical velocity. The vortex and antivor- 
tex are created at the lateral sides of the disk, and they 
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FIG. 1: Dynamics of the density 



2 /no and phase 



= arg^ profiles for the Gaussian potential with Vo/(gno) = 
—5, v/v s = 0.6, and d/£ = 4.5 in the frame comoving with the 
potential. The arrows indicate the directions of circulations. 
The field of view of each panel is 160£ x 67£. 



are released into the wake (see Fig. 1 of Ref. The 
Gaussian potential dU) for Vo > also exhibits similar 
behavior 

By contrast, we investigate the dynamics for a red- 
detuned (Vo < 0) moving Gaussian beam, which is shown 
in Fig. [TJ First, a density dip appears at the head of the 
potential [Fig.[T](a)], which splits into vortex and antivor- 
tex [Fig. Q] (b)]. The vortex pair then passes through the 
potential avoiding the high density region [Fig. [1] (c)], 
and it is released behind the potential. The released 
vortex dipole moves in the — x direction with the ve- 
locity h /(ml) = v 3 £/l, where I is the distance between 
the vortex and antivortex. Such vortex dynamics are re- 
peated and vortex dipoles are shed in the wake periodi- 
cally. The train of vortex dipoles is dynamically unstable 
and eventually sinuous destabilization spreads it in the 
±y directions [13, [HI • In contrast to the case of a blue- 
detuned Gaussian beam, the Benard-von Karman vortex 
street [17| was not found for a red-detuned beam to the 
best of our investigation. 

For an attractive potential (Vo < 0), a vortex dipole is 
nucleated at the head of a moving potential, whereas a 
vortex dipole is nucleated at the lateral sides for a repul- 
sive potential (Vo > 0). To understand these behaviors, 
we consider the irrotational flow of an inviscid and in- 
compressible fluid around a rigid disk. The velocity field 
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FIG. 2: Critical velocities for vortex dipole creation versus 
beam waist d of a moving red-detuned Gaussian beam. When 
the velocity of the beam is adiabatically increased, the vor- 
tex dipole creation occurs at v C ri (circles). When the velocity 
of the beam is adiabatically decreased, the vortex dipole cre- 
ation stops at v C r2 (squares). The intensity of the Gaussian 
potential Q is fixed to Vo/(gno) = — 5. 



in the frame comoving with the disk is given by [2c 

R 2 _o 



U 1 



(5) 



where U and R are the velocity and radius of the disk, 



(x 2 



y 



2 ) x / 2 , and — arg(x + iy). The velocity 



becomes maximum at (r, 0) = (R, ±7r/2), which agree 
with the locations at which vortices are nucleated. This is 
because vortices are nucleated when the velocity exceeds 
the local Landau critical velocity. Moreover, the density 
is low where the velocity is large due to Bernoulli's law 
and a quantized vortex with a density hole is easy to be 
nucleated. 

For a disk-shaped strong attractive potential, where 
the density inside the disk n- in is much larger than the 
outside density n out , the inside velocity must vanish in 
the limit of ni n /n out —> oo since the normal compo- 
nent of the flow must be continuous across the boundary, 



J 



fWt^cmf Therefore, the tangential component 
of the velocity must vanish at the boundary, since it 
must be continuous. Such velocity field outside the disk 
is given by 



R 2 

iv y = U ( H ^e" 



(6) 



The velocity becomes maximum at r = R and = 0, 7r, 
and thus a vortex dipole is nucleated at the head of a 
potential. Perturbations at the tail of a potential (r = R 
and = 0) do not grow because they go downstream 
away. 

It is interesting to note that the system exhibits hys- 
teresis behavior with respect to the velocity. When the 
velocity of the beam is adiabatically increased from zero, 
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FIG. 3: (a)-(h) Snapshots of the density \ip\ 2 /no and phase 
<j) — arg^ profiles in the rest frame. A Gaussian potential 
with Vo/(gno) = —5 is moved as in Eq. J7} with vo/v s = 0.47 
and T It — 17. The beam waist is d/^ = 3.54 in (a)-(d) and 
d = 2.67 in (e)-(h). The arrows in (c) and (g) indicate the 
traveling direction of the vortex dipoles. The field of view of 
each panel is 78£ x 33£. (i) Velocity of the launched vortex 
dipole versus beam waist d. The parameters except d are the 
same as those in (a)-(h). 



the successive vortex dipole creation starts at a critical 
velocity v cr i, which is shown by the circles in Fig. [2j On 
the other hand, when the velocity of the beam is de- 
creased from v > v cr i, the vortex dipole generation does 
not stop at v = v cv \. The velocity v cr 2 (squares in Fig. [2]) 
at which the vortex dipole generation stops is slower than 
v cr i . For the velocity v cr 2 < v < v cr i , the stationary flow 
remains to be stationary, and the successive vortex dipole 
creation is kept once it is triggered by perturbations. 



B. Launching a vortex dipole 

Using a red-detuned Gaussian laser beam, we can 
launch vortex dipoles in a controlled manner. The initial 
state is the ground state for the Gaussian potential in 
Eq. (j4j) with Vo < and v = 0. The potential then starts 
to move in the — x direction at t = and stops at t = T 



4 



as 



(t<0,t>T) 
v (0 < t < T) 



(7) 



Figures [3] (a)-|3] (h) show launching dynamics of a vortex 
dipole in the rest frame. After the stop of the potential, 
the created vortex dipole is released from the potential 
and continues to move in the — x direction. Figures [3] 
(e)-|3] (h) show the case of a smaller potential, in which 
a vortex dipole is bound so tight that the topological 
defects cannot be discerned [Fig. 03 (h)]. Such a solitonic 
structure is referred to as a rarefaction pulse [27j, which 
has been shown to be formed by a local depletion [28| 
and an oscillating potential [2l[. The rarefaction pulse 
propagates faster than the vortex dipole [Figs. [3] (c) and[3] 
(g)]. Figure[3](i) shows the dependence of the velocity v v a 
of the launched vortex dipole on the beam waist d. Since 
the distance between the created vortex and antivortex 
increases with an increase in the size of the potential, the 
velocity of the vortex dipole is a decreasing function of 
the size of the potential. Thus, the vortex dipole can be 
launched with controlled velocity. 



IV. VORTEX GENERATION BY MULTIPLE 
GAUSSIAN BEAMS 

A. Flow around double beams 

We next consider the case of multiple Gaussian po- 
tentials. For simplicity, we assume that two identical 
Gaussian potentials move at the same velocity as 
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(x-Xj+vt) 2 + (y-yjY 
d 2 



(8) 



The positions of the potentials are, without loss of gen- 
erality, (xj,yj) — (T§ si n X? if cos x) with the upper 
signs for j = 1 and the lower signs for j = 2, where the 
angle x determines the alignment of the potentials. 

Figure H] shows a typical dynamics for two blue-detuned 
(Vq > 0) beams with % = 0. The initial state is the 
ground state for v = and the potentials start to move 
in the — x direction at t = 0. After the initial disturbance, 
the system exhibits periodic behavior. A vortex dipole is 
created between the potentials and released in the wake 
[circle in the left panel of Fig. [4] (a)] with its counter part 
being remained in the potentials [circles in the right panel 
of Fig. |H(a)]. Then, the topological defects remained in 
the potentials are released from both ends [Fig. 2] (c)], 
and this dynamics is repeated. We note that small vortex 
dipoles released from between the potentials propagate 
in the x direction, whereas vortices in the outer rows 
propagate in the — x direction [arrows in Fig. 01(e)] . 

In order to understand the behavior in Fig. 3J we con- 
sider the irrotational flow of an inviscid and incompress- 
ible fluid around two rigid disks with radius R > D/2 
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FIG. 4: Dynamics of the density \i/j\ 2 — \ip\ 2 /no and phase 
4> — arg^ profiles for two Gaussian potentials with Vb / (gno) = 
100, d/i = 1.58, v/v s = 0.35, and D/£ = 12.6 in the frame 
comoving with the potential. The circle in the left panel of (a) 
indicates a vortex dipole released from between the potentials, 
and those in the right panel its counterpart remained in the 
potentials. The arrows in (e) indicate the directions in which 
the vortices propagate. The field of view of each panel is 
160£ x 67£. 



located at (x,y) = (T^i snl X? ±y cos x)- The velocity 
field is given by [2J-[2S| 



(z/i?) 2 + sinh 2 /3L 



e %x p(ir - if}) 



+e lx p(ir + if}) + —CM cos X 

7T 



(9) 



where f} = cosh 1 [D/(2R)], z = x + iy, r = f3 + 
\og(ie~ %x z/ R + sinh f}) — \og(ie~ %x z / R — sinh/3), and p 
and C are Weierstrass's elliptic functions with half peri- 
ods 7T and i cosh" 1 [(D 2 /R 2 - 2) /2]. The velocity field ® 
approaches to v x — iv y = U for \z\ — » oo. For x — 0? 
which corresponds to the case in Fig. HJ the velocity field 
in Eq. (|9j) is maximum at the inner edges of the disks. We 
denote the maximum velocity as v max = v x (x = 0,y = 
±D/2 =p R). The contour plot in Fig. [5] shows v max as 
a function of the distance D between the disks and the 
unperturbed velocity U. The maximum velocity v max 
monotonically decreases with an increase in D and lin- 
early increases with an increase in U. 

In Fig. [5j we superimpose the plots of the critical ve- 
locity v cr for vortex dipole creation, obtained from the 
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FIG. 5: Contour lines show the maximum velocity v ma x as 
a function of D/R and U/v a , obtained from Eq. ((9]). The 
contour lines indicate v m ax/v s = 0.05, 0.1, 0.15, 0.2, • • • from 
below. The plots show the critical velocity v cr for vortex 
dipole nucleation versus D / i? e ff , which is obtained by simula- 
tions as in Fig. [U Here i? e ff is the radius of the density hole 
produced by the potential. The beam waist and the effective 
radius are Reft/£) — (5,9.5) for the circles and (10,18) 
for the squares. 



GP equation (|3j) numerically. In the time evolution, we 
adiabatically increase the velocity v of the potential from 
v = until a vortex dipole is nucleated as in Fig.HJ which 
determines the critical velocity v cr . From Fig. [5j we find 
that the plots roughly fit the contour lines, implying that 
a vortex dipole is nucleated when the maximum velocity 
between the disks exceeds a local critical velocity, which 
is almost independent of D. 



B. Launching a vortex dipole 

Using two blue-detuned Gaussian beams, we can 
launch a vortex dipole, as in the case of a red-detuned 
beam in Sec. lIIIBl The initial state is the ground state in 
the presence of two Gaussian potentials (|5J) with Vb > 0, 
v = 0, and x = 0- The potentials are then displaced as in 
Eq. ([7]). Figures |5](a)-|5|(c) show the launching dynamics 
in the rest frame. A vortex dipole is generated between 
the potential and moves in the direction opposite to that 
of the displacement of the potentials. A vortex dipole is 
thus launched by "pulling" the potential, which is in con- 
trast to the case of a red-detuned beam, in which a vortex 
dipole is launched by "pushing" the potential. After the 
launch, a counterpart of the released vortex-antivortex 
pair remains at the potential [dashed circle in Fig.[6](d)], 
which prevents the successive launch. By contrast, the 
successive launch is possible for a red-detuned beam. 

Figure [6] (e) shows the dependence of the velocity v v( ± 
of a launched vortex dipole on the distance D between 
the two Gaussian potentials. The velocity v v d decreases 
with an increase in D, since the distance between the 




FIG. 6: (a)-(c) Dynamics of the density profiles \i/j\ 2 /no for 
two Gaussian potentials with Vo/(gno) = 5, d/£ = 4.5, and 
D/£ — 19 in the rest frame. The potentials are moved as in 
Eq. © with vo/v s = 0.24 and T/r = 133. The arrows in 
(a) and (c) indicate the direction of the displacement of the 
potentials and that of the propagation of the vortex dipole, 
respectively, (d) Phase profile, arg^, at t = t/r = 220. The 
topological defects in the dashed circle are the counterpart of 
the launched vortex dipole. The field of view of each panel 
is 90£ x 65£. (e) Velocity of the launched vortex dipole ver- 
sus distance D between the two Gaussian potentials. The 
parameters except D are the same as those in (a)-(d). 



launched vortex and antivortex increases. 



C. Trapping of a vortex dipole 

We consider the dynamics for two red-detuned beams 
as shown in Fig [TJ The initial state is the ground state 
for two Gaussian potentials in Eq. J5J) with Vb < 0, v = 0, 
and x — [Fig [3 (a)] and the potentials start to move in 
the —x direction at t = 0. Initially, the sudden increase 
in the velocity disturbs the system and several vortex 
dipoles are generated near the potential [Fig 13(b)] . The 
vortex generation is only at the starting, and the gen- 
erated vortex dipoles depart from the potential except 
the one between the potentials. Interestingly, this vortex 
dipole is trapped between the potentials and oscillates 



(e) t=5.4 




o 



FIG. 7: Dynamics of the density \tp\ 2 = \^\ 2 /no and phase 
4> — arg^ profiles for the potential in Eq. ((8} with Vo/(gno) = 
—5, v/v s = 0.6, <Z/£ = 4.5, and D/t; = 31.6 in the frame 
comoving with the potential. The field of view of each panel 
is 160£ x 67£. 



around the center in the moving frame [Figs [7] (c)-|7] (e)]. 

In order to understand the trapping of a vortex dipole 
between moving potentials as shown in Fig.[7J we perform 
linear stability analysis assuming inviscid, incompress- 
ible, and irrotational flow. We consider the problem in 
the frame comoving with the potentials. Two point vor- 
tices with circulations —h/m = —tz and k are located at 

Oi,?/i) = (£1,2/2 + 771) and 02,2/2) = (6,-V 2 + ^), 

where £1, 771, £2, an d 772 are small deviations from the 
stationary state and I is the distance between the vor- 
tex and antivortex for the stationary state. Equations of 
motion for the deviations are given by (j = 1,2) 



S,j — Vx {ptj 1 Uj ) 



2tt r 2 ' 

K X\ — X2 
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(10a) 
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where r 2 = (x x - x 2 ) 2 + (yi - y 2 ) 2 = (£1 - &) 2 + (I + ^1 - 
r]2 ) 2 . The first terms on the right-hand side of Eq. (fT0|) 
are the velocity field without the point vortices and the 
second terms are the velocity field produced by the other 
point vortex. We impose the same boundary condition 
as in Eq. (|6j) at x 2 + (y±D/2) 2 = R 2 , giving the velocity 



FIG. 8: Stability of a trapped vortex dipole with respect 
to the distance D between the potentials and the normalized 
velocity U = 2tv RU/k. When A in Eq. ([T3]) is real, the system 
is dynamically unstable. The contour plot shows 2tyR 2 \X\/k 
where A is pure imaginary. The circle corresponds to the case 
in Fig. [3 



field without the point vortices as 
4/7sinh 2 /3 



{z/Rf 



-—2— [p(ir - i/3) - p(ir + i/3)] , 
smh p 

(11) 



where /3 and r are defined in Sec. IIV A[ Substituting 
Eq. ([TT]) into Eq. ([TO]) and neglecting the second and 
higher orders of the small deviations, we have 
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(12b) 
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where = v x (0,±l/2) and c = —d y v x (0,l/2) = 
d y v x (0,-l/2) = -d x v y (0,l/2) = d x v y (0,-l/2). In order 
for the point vortices to be stationary at (0,±Z/2), we 
set vo = «/(27rZ), which determines Z for given J7 ? D, and 
i?. Equations (fT2j) then become homogeneous differential 
equations and have the solutions of the form Cj^Vj ^ eXt - 
The eigenvalues A are obtained as 



A = ± W c 2 -r 



rt 2 ' 



(13) 



If the eigenvalues are pure imaginary, i.e., c < n/{i:l 2 ) 
the solutions are oscillatory and the two point vortices 
are stable. If c > ft/(7r/ 2 ), the small deviations grow 
exponentially and the system is dynamically unstable. 

Figure [8] shows a contour plot of |A| in Eq. (jT3j) for the 
stable region. The stable trapping dynamics shown in 
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FIG. 9: Stroboscopic images of the density profile \ip\ 2 /no 
in the rest frame. A vortex dipole is launched by the two 
Gaussian potentials as in Fig. [6] with Vo/(gno) = 5, d/£ — 
4.5, and D/£ = 19. The stroboscopic images are taken at 
time intervals of At/r = 20. (a) A Gaussian potential with 
Vo/(gno) = 1 and d/£ = 4.5 is located at y = 1.6. (b) A 
Gaussian potential with Vo/(gno) = —1 and d/£ = 4.5 is 
located at y — 1.6. (c) A Gaussian potential with Vo/(gno) = 
1 and d/£ = 4.5 is switched on at t/r = 330. (d) A Gaussian 
potential with Vo/(gno) — —1 and d/£ = 4.5 is switched on at 
t/r = 330. The field of view of each panel is 224£ x 80£. (e), 
(f) Schematic illustrations for explaining the Magnus effect. 



Fig. [71 corresponds to the circle in Fig. [51 From Fig. [5J 
we find that |A| increases with an increase in U and the 
trapping is tight for large U. However, if U exceeds the 
critical velocity, vortex dipoles are successively generated 
from the potentials as in Fig. [TJ which disturbs the trap- 
ping. Therefore, there is an appropriate region of D and 
U for the most efficient trapping. 



V. MANIPULATION OF VORTEX DIPOLES 

We next show that the velocity and trajectory of a 
launched vortex dipole can be changed after the launch. 
We use a static Gaussian potential to curve trajectories 
of vortex dipoles, and use a time-dependent Gaussian po- 
tential to accelerate and decelerate them. Figure [9] shows 
stroboscopic images of such dynamics. A vortex dipole is 
launched in the same manner as in Fig. [6] In Fig. [9] (a), 
a weak repulsive Gaussian potential is located near the 
trajectory of the vortex dipole and the trajectory bends 
toward the potential. In Fig. [9] (b), the obstacle poten- 
tial is attractive and the trajectory bends in the opposite 
direction. In Fig. [9](c), a repulsive Gaussian potential is 
switched on at the moment that the vortex dipole passes 
through, which accelerates the vortex dipole. In Fig. [9] 
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FIG. 10: Collisions between vortex dipoles launched by the 
method in Fig. [3] The directions in which the vortex dipoles 
propagate are indicated by the arrows. The numbers and 
their colors identify the vortices, (a) Collision with the impact 
parameter b/£ = 3.16. Vortex exchange occurs, (b) Collision 
with b/£ = 6.32. The two vortex dipoles pass each other, (c) 
Collision of vortex dipoles with different velocities. The beam 
waist is d/£ = 3.54 for the left-hand side in (c) and d/£ = 3.16 
for the others. The other parameters are the same as those 
in Fig. [3] The field of view of each panel is 78£ x 33£. 

(d), on the other hand, an attractive potential is switched 
on and the vortex dipole is decelerated. 

The behaviors in Figs. [9] (a)-|9] (d) can be understood 
by the Magnus effect. The Magnus force on each vortex 
is given by Fm = mnn x Av, where k is the vector of 
circulation and Av is the velocity relative to the unper- 
turbed velocity of the vortex dipole. When force F is 
exerted on a vortex, the velocity changes in such a way 
that F balances with Fm, and hence 

Av x F. (14) 

In Fig. [9j the origin of the force is the density gradient 
produced by the obstacle Gaussian potential. The en- 
ergy of a vortex is roughly given by E ~ 7rh 2 n/m and 
the force on a vortex is F ~ —irk 2 Vn/m, which is in the 
direction opposite to the density gradient. Two exam- 
ples of the effect of density gradient on a vortex dipole 
are shown in Figs. [9] (e) and [9] (f). When the density 
gradient is perpendicular to the traveling direction of a 
vortex dipole [rightward for Fig. [9] (e)], its trajectory is 
curved toward the direction of — Vn, which explains the 
dynamics in Figs. [9] (a) and[9](b). When a vortex dipole 
propagates in the direction of Vn as in Fig. [9] (f), the 
distance between the two vortices decreases and the vor- 
tex dipole accelerates. This corresponds to the situation 
in Fig. [9] (c), where a density dip suddenly appears just 
behind the vortex dipole. The deceleration in Fig. [9] (d) 
can also be understood in a similar manner. 

Using the method to launch vortex dipoles, we demon- 
strate collisions between two counter-propagating vortex 
dipoles. The two vortex dipoles are launched in the ±x 
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directions by red-detuned Gaussian beams as in Fig. [3l 
Figures [10] (a) and [10] (b) show collisions between two 
vortex dipoles launched with the same velocity. When 
the "impact parameter" b is small, each vortex in the 
two vortex dipoles are exchanged, and the two vortex 
dipoles with new pairs go away with the same scattering 
angle [Fig. [10] (a)]. The scattering angle is tt/2 for frontal 
collision (b = 0, data not shown). In Fig. ITUl (b). the im- 
pact parameter is large and the two vortex dipoles pass 
each other without vortex exchange. Figure [10] (c) shows 
collision between vortex dipoles with different velocities, 
where the vortex exchange occurs and the new vortex 
dipoles scatter with scattering angles ip and n — <p. This 
is understood from the momentum conservation, where 
a slower vortex dipole (larger distance between vortices) 
has a larger momentum [18[. Collisions between point 
vortex di pol es in classical fluids have been studied in 
Refs. [H, Hjl . 

VI. CONCLUSIONS 

We have investigated dynamics of a 2D BEC stirred 
by far-off-resonant Gaussian laser beams. When a red- 
detuned beam (attractive potential) moves through a 
BEC exceeding a critical velocity, a vortex dipole is nucle- 
ated at the head of the moving potential (Fig. [TJ. This is 
in contrast to a blue-detuned (repulsive) beam, in which 
a vortex dipole is generated on both sides of the poten- 
tial. The vortex shedding from a red-detuned beam ex- 
hibits the hysteresis behavior with respect to the velocity 
(Fig. [2]). When a red-detuned beam is displaced, a cre- 
ated vortex dipole propagates towards the direction of 
the displacement, which can be used as a "vortex-dipole 



launcher" (Fig. [3]). When two repulsive potentials are 
moved, a vortex dipole is generated between the poten- 
tials, which propagates to the direction opposite to the 
motion of the potential (Fig. [4]). A vortex dipole can 
also be launched by two blue-detuned beams (Fig. [6]). 
The velocity of vortex dipoles can be controlled in the 
two launching methods. For two red-detuned beams, we 
found that a vortex dipole can be trapped between two 
moving potentials (Fig. [7]). Launched vortex dipoles can 
be manipulated by an external potential (Fig. EJ). A tra- 
jectory of a vortex dipole is curved by a Gaussian po- 
tential located near the trajectory. A Gaussian potential 
with time-dependent intensity on the trajectory of a vor- 
tex dipole can change its velocity. 

We have thus shown that we can create and launch vor- 
tex dipoles in a controlled manner and also manipulate 
their trajectories and velocities after the launch. These 
are achieved only by a quasi-2D BEC and Gaussian laser 
beams. This technique enables us to investigate a variety 
of dynamics of quantized vortices in superfluids, such as 
collisions of vortex dipoles (Fig. [TO]). 
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